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The semileptonic decays r~ — > fi~M(M — tv°, r\, r\ ) could be sensitive probe for new physics sce- 
narios with lepton flavor violation(LFV). Motivated by the recent Belle measurement, we investigate 
these decays in type III two-Higgs-doublet model (2HDM III), R-parity violating supersymmetric 
\Q ' models(RPV SUSY) and flavor changing Z' models with family non-universal couplings, respec- 

o ' 

tively. In these new physics scenarios, there are LFV couplings at tree level. Our results have 
^ , shown that the decays are very sensitive to the LFV couplings and could be enhanced to the present 

Fj^ ' experimental sensitivities. We have derived strong bounds on relevant couplings of these models, 

which may be useful for further relevant studies. 
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I. INTRODUCTION 



The flavor physics of fermions are among the most mysterious fundamental problems in particle physics. In the 



Qh. standard model(SM), neutrinos are exactly massless due to the absence of the right handed chiral states (vr) and 



X 

H 

observation 



isci 

.1 



the requirement of SU(2)l gauge invariance and renormalizability, so that the chirality conservation implies lepton 
number conservation. In the past decade, the most exciting progress in understanding of these issues has been the 
of oscillation of atmospheric neutrinos with very large mixing. The observation shows that neutrinos 
are massive and the lepton flavor violating (LFV) exists in the neutral lepton sector. 

In the SM supplemented with massive neutrinos, the neutrino mixing will induce, at loop level, rare LFV processes 
in charge lepton sector such as r — > /17, r — > /xM(M=light hadrons), etc. These processes are expected to be 
proportional to the ratio of masses of neutrinos over the masses of the W bosons, which is negligible small. However, 
the t — » fx mixing could be large in new physics models 2], thus the LFV r decays provide some sensitive probes for 
new physics beyond the Standard Model. 

It is interesting to note that the B factories BaBar and Belle arc also t factories. The r production cross section at 
BaBar and Belle is as large as cr e + e -_^ r - T + = 0.89nb. The integrated luminosity up to now at Belle is about 540/6 -1 , 
which corresponding to about 4 x 10 s t~t + pairs. While integrated luminosity at BaBar is about 320/6~ 1 . As BaBar 
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and Belle are steadily accumulating more data, it would be very promising to search for rare r decays at BaBar and 
Belle to constrain or reveal new physics effects. 

Ex „ ^ ^ bee „ Peltad ,o, , „ re d _ _ LFV ^ - , eM , flUQ, 
0, r -> e(/x)7r+7r- Q 



and r — > e(fi)j 



The primary theoretical studies are focused on zi and r radiative 

nn 

decays and their decays to three charged leptons in different models 9j, llOj . There are also many studies on decays 
t — > IM. llakovac et al., have studied r — * e(fi)M decays in models with heavy Dirac or Majorana neutrinos and 
found t — * e(f>(p Q , n°) with branching ratios of order f CP 6 1 lj . Sher has investigated the decays r — ► fin in the large 



tan region in seesaw MSSM 

r~l 
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In MSSM model, these decays are analyzed by Fukuyama et al [12j. Brignolc 
and Rossi [13 has presented a comprehensive theoretical study on these decays recently in a general unconstrained 
MSSM. A model independent study has been preformed by Black et al [lj] to bound new physics scales. 

In order to get information on LFV couplings, an important observable, the anomalous magnetic moment of the 
muon [g— 2) M , should be considered additionally. Several works constraining LFV processes from previous estimations 
on a M have been carried out in Refs|l J. 17|. 

The most recent experimental search for semileptonic LFV r transitions has performed by Belle [4| using only 
153.8/6- 1 data 

B(t -> /i7T°) < 4.1 x 10~ 7 , 90% CL 
B{t fir]) < 1.5 x 1(T 7 , 90% CL 

B(t -> fir)') < 4.7 x 10~ 7 , 90% CL. (1) 



These results are already 10 to 64 times more restrictive than previous CLEO limits|3|]. Using the theoretical formula 
derived by Sher u| in a seesaw MSSM, the Belle new results have improved the constraints of the allowed parameter 
space for rriA — tan fj in MSSM. It would be very worthy to study these decays in other new physics models to derive 
bounds on relevant parameters. In this paper, we will study these decays in three new physics scenarios, namely, 



• the 2HDM III 18, 19, 20], where FCNC and LFV could arise at the tree level in the Yukawa sector when the 
up- type quarks and down quarks are allowed simultaneously to couple more than one scalar doublet, 



n Ei. 3 3 



• SUSY theories with R-parity broken |2lL 122 . 23] , in which the R parity odd interactions can violate the lepton 
and baryon numbers as well as couple the different generations or flavors of leptons and quarks, 



• flavor changing Z 1 models 



E3 3 



27 1 with family non- universal couplings. 



We have shown that the LFV semileptonic r decays could be enhanced to the present B factories sensitives in 
the above three scenarios, thus we have derived bounds on the LFV couplings in the models, which are tighter than 
exiting ones in the literature. 



3 

In next section, we present calculations of the decays and bounds on the LFV couplings in the aforementioned three 
scenarios. Finally in Sec. Ill, we give our conclusions. 



II. MODEL CALCULATIONS 



A. Hadronic matrices of local operators 



Before detail model calculations of r — > \xM(M = 7r°, r\, rf), we would specify hadronic matrices elements which are 
inputs for calculating these decay amplitudes. 
At first, we need 



(n°(p)\q^ 75 q\0)= -j=UPm 
with = 130 ± 5 MeV, and the so-called chiral condensation matrix 



(2) 



(3) 



where q=u or d. We will see that the l/2m q factor will cancel the corresponding quark mass coupling of the weak 
scalar interaction operators, and thus enhances scalar interaction contributions. 

As for 77 and 7/, the situation is much more complicated than it . The relevant matrices are defined by 



{M(p)\s^ 5 s\< 



V2 



2m q (M(p)\q 75 q\0) 



— h q 

V2 M ' 



2m s (M(/j)|s75s|0) = -ih s M , 



(4) 



where M — rj or 7/. These equations define eight non-perturbative parameters, which however are not all independent. 
They could be related to fewer independent non-perturbative par ameters by 77 — rj' mixing scheme. In this paper, 
we will take the Feldmann-Kroll-Stech (FKS) mixing scheme 30] . In FKS mixing scheme the parton Fock state 
decomposition can be expressed as 



/ 



I''/) 



/ 



V \rf. 



sm ( 



sin m cos ( 



V \Vs) 



(5) 



where <j> is the mixing angle, \r) q ) — (\uu) + \dd)) /y/2 and |ry s ) = \ss). The four parameters P M are therefore related 
by 



/« = f q cos </>, /" = -f s sin < 



= fq Sm <t>> /-)' = C0S <f>> 



(6) 



and an analogous set of equations for the h 



M 



hi = h q cos <ft, h^ = —h s sin < 



Moreover the ft g and h s could be related to f q , f s and mixing ang le <j> Q 

ft 9 = iq (m 2 cos 2 + m 2 / sin 2 0) — (m 2 , — m 2 ) sin </> cos </>, 

/i s = / s (ml, cos 2 + m 2 sin 2 cf>) — ^j=(m 2 , — m 2 ) sin</>cos0. (8) 

V2 

The three remaining parameters / 9 , / s and 4> m FKS scheme have been constrained from the available experimental 
data with results 30j 

f q = (1.07 ± 0.02) /„, f s = (1.34 ±0.06)/,,, cf> = 39.3° ± 1.0°. (9) 

From these parameters and /„•, Beneke and Neubertjiu have derived 

/« = 108 ± 3MeV, f° = -111 ± 6MeV, ft« = 0.001 ± 0.003Gel/ 3 , ft* = -0.055 ± 0.003GeV 3 , 

= 89 ± 3MeV, = 136 ± 6MeV, h\, = 0.001 ± 0.002G e y 3 , ft* = 0.068 ± 0.005Gel/ 3 . (10) 

It should noted that ft^ and ft^, are poorly determined. In our numerical calculations, we take these hadronic 
parameters with la variant to display their uncertainty effects on our bounding on LFV couplings. Now we are ready 
to calculate the decays in the aforementioned three new physics scenarios. 



B. r -> fiM decays in 2HDM III 

Two Higgs Doublet Model(2HDM) is the popular and the most simplest extension of the SM with a scalar sector 
made of two instead of one complex scalar doublets. In order to build a 2HDM without FCNC at tree level, it is 
achieved by requiring either that w-type and c?-type quarks couple to the same doublet (2HDMI) or that u-type quarks 
couple to one scalar doublet and d-type quarks to the other(2HDMII). The case in which scalar FCNC not forbidden 
is dubbed 2HDMIII, where the Higgs doublets could couple to both the u- and d type quarks at the same time 2, 2p| . 

Generally one can write Yukawa Lagrangian of 2HDMIII QflD 



Cy = vYA.lHiUj.r + '/.'Qu.ll.lh.u + ^/Ij.HJ,.,, + ^Qi,LH 2 D jtR + h.c, (11) 

where Hiii = 1,2) are the two Higgs doublets. Qi h is the left-handed fermion doublet, Uj } n and DjR are the 
right-handed singlets, respectively. These Qi t L, Uj^R and Dj.R are weak eigenstates, which can be rotated into mass 
eigenstates, while ri U ' D and £ ' are the non-diagonal matrices of the Yukawa couplings. 

For convenience one can express Hi and H2 in a suitable basis such that only the r^- D couplings generate the 
fermion masses, i.e., 

I \ 

m = 0. (12) 




The two doublets in this basis have the form 



Ht = 



1 

71 



o 



Hi 



1 

71 



<#j + iA° 



(13) 



72 G+ 

where G 0,± are the Goldstone bosons, H^ and A are the physical charged-Higgs boson and CP-odd neutral Higgs 
boson respectively. The advantage of choosing the basis is the first doublet H\ corresponding to the scalar doublet 
of the SM while the new Higgs fields arising from the second doublet H 2 . So H% does not couple to gauge bosons of 
the form H%ZZ and H%W + W~ . 

In Ea. l|13[) . 4>i an d ^2 are n °t the neutral mass eigenstates but linear combinations of the CP-even neutral Higgs 
boson mass eigenstates H° and h° 



H° 



h° = 



sm a 



, cos a, 



(14) 
(15) 



where a is the mixing angle. In the case of a = 0, (0°, </>§) coincide with the mass eigenstates of H° and h°. 
After diagonalizing the mass matrix of the fermion fields, the Yukawa Lagrangian becomes 

Cyc = ///:'■',,/,■ + vFjQi,LHiD jjR + i,',^, / //,',/,' + !; ! /,(>,j.I I ■!),.!< + h.c, 
where Qi,L, Uj^n, and Dj r now denote the fermion mass eigenstates and 



nrU,D\-l U,D j r U,D ■ D r U,D\ 

( v l ) -V ' -V R — —M {Mi/ =S ljmj ' ), 



(16) 

(17) 
(18) 



In Eq.©, V£ are the rotation matrices acting on up and down-type quarks, with left and right chiralities respec- 
tively. Thus Vckm = {YlV^l ^ s the usual Cabibbo-Kobayashi-Maskawa (CKM) matrix. We can see from Ea. llGI) 
that the matrices £ U ' D , as defined by Ea. 1)18(1 . allow scalar-mediated FCNC. That is, in the quark mass basis only 
the matrices f] U ' D of Ea. (|17fl are diagonal, but the matrices £ U ' D are in general not diagonal. The FCNC part of the 
Yukawa Lagrangian is 



C 



H° sin a + h° cos a 



Y.FCNC 



V2 



U 



^(1 + 75) +^^(1-75) 



U + D 



e D i(i + 75) + e Dt i(i-75) 



D 



71 



U 



^^(1 + 75) -^^(1-75) 



U-D 



l D ^(i+75)-e Dt ^(i-7 5 ) 



3 3 



D 



(19) 



The corresponding Feynman rules from Ea. (|19fl can be found in Refs. 

Because the definition of & D couplings is arbitrary, we can take the rotated couplings as the original ones and 
shall write £ U,D in stead of £ U,D hereafter. 



tU,D _ \ . 

Si 7 — A n 



(20) 



which ensures that the FCNC within the first two generations are naturally suppressed by small fermions masses. This 
ansatz suggests that LFV couplings involving the electron are naturally suppressed, while LFV transition involving 
muon and tau are much less suppressed and may generate sizeable effects. 
The decay amplitudes are given by 



A(t- -> /i-7r°) 
A{t~ -> iTrj) 
A(t- 



G f 

— y/m^mrUomlo {EX t ^(Jit) s+ p + FX*^(p,T) S -p} , 

Gf 
2~7I 

Gf 
2V2 



^fm~m~ T 



h\J q + h s v ,3 s 



Kh{v>t)s+p + 



h%K q + h°K a 



h\K q + h s v ,K s 



where h % M are defined by Eq0| and the auxiliary functions are 

1 / , \ ./ sin 2 a cos 2 a\ , 

E = I ReXuu + ReXdd +i I — s— + — — 1 I M„ 



F = — 



"A 
1 



rn 



H" 



it) 



ImXdd I , 



m 



sin 2 a cos 2 (x 

2— . ReX uu + ReXdd ) +i [ — 5 1 5 — ) I Im ^uu - ImX d d 



J q = 



J, = 



K q = 



A° 
1 



m 



H" 



m 



h° 



m u + m d 
1 



— i— (m u ReX uu - rridReXdd) +i ( - + C ° S ? - ) ( m u ImX llu + mdlmXdd 



-ReX s 



sin 2 a cos 2 a 
«l — s 1 s — JmA» 



'7i° 



1 



K, 



m u + m,d 
1 



A" 



m u ReX uu - m d ReXdd +i 



• 2 9 
sm a cos a 



'/!» 



m u ImX uu + rridlmX 



<i,i 



-ReX ss + i 



m 



A 



■ 2 o \ 

sm a cos a 



ImXs 



(21) 
(22) 
(23) 



(24) 



From Ea 12 11231 we can see that there are two types operators (/j,t)s±p contribute to decays at tree-level. And it 
should be noted that the small quark mass factor in the neutral Higgs- quark couplings are cancelled by the one in 



TABLE I: Constraints on the A rfl from r — > ^i7r° , firf- ' decays 



Decay modes Bounds on X Tf , 



Previous Bounds 



T — » fnv 

T -> (Iff 



< 63.62 

< 2.76 

< 4.48 



X T fj, ~ 0(1) [2,8], A TU ~ 0(10) - O(10 2 ) [24J 
A TM ~ O(10 2 ) -o(io 3 ) [25] 
|A MM | = A TT = I A Mr = |A eM | = 10 [36] 



In this work we choose £ ' D to be complex for the sake of simplicity, so that besides Higgs boson masses, only 
A MU , Add and X ss in the quark sector and A ru in the lepton sector are parameters relevant to the semileptonic r decays. 
Totally there are seven parameters in the amplitudes: X qq (q = u, d, s), A ru , their phases 9, the masses of neutral Higgs 
myfi , mjio , m^o and the mixing angle a. Taking into account of constraints on parameters from experimental data 
and theoretical limits, the values of parameters can be taken as 



)i constraint 



6 = tt/4, |A mi | = 150, |A d d| = 120, |A„| = 100, m h o = H5GeV, m A o = UOGeV, m H o = WOGeV, a = tt/4.(25) 
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FIG. 1: The branching ratios as functions of A Tfl for fixed m H o = 150GeV, m A o = 120GeV,m h o = 115Gev,a — 7r/4, the 
phase of A TM , 8 = tt/4 in the 2HDM III. (a) for r~ —> /i~7r° decay, (b) for r~ — > /^"r/ decay, and (c) for r~ — > /U _ r?' decay. The 
upper horizontal lines denote Belle upper bounds. The lower horizontal lines denote the possible upper bounds obtained naively 
by scaling Belle upper bounds with 540/6 -1 integrated luminosity. The dotted bands correspond to the effects of hadronic 
inputs uncertainties. 



where the higgs masses satisfy the relation WbGeV < m h o < m^o < m H o < 200GeV 

In the literature , th e effects of Ay have been studied under different phenomenological considerations 



3 3J3 



3jj|37j]. For comparison we list bounds on the X TfJ _ parameter of 2HDM III in Tabled. In Rcf. 



y, 



y, 



it is suggested X Tfl ~ 0(1). However, based on limits from (g — 2) M results y4j and tau decay [25j, X TfJ , could have 
much larger value ~ 0(10 — 10 3 ). The branching ratios for t~ — > fi~ir°(ri, rf) versus \X Tfi \ in 2HDM III are presented 
in Fig.l, where Fig.l.(a), Fig.l.(b) and Fig.l.(c) are the results for t~ — > /i~7r°, r~ — > yTr) and r~ — > /i~r/ decays, 
respectively. The dotted bands show our theoretical uncertainties due to the variants of hadronic inputs as in Eci. H10[) . 
The hadronic uncertainties in B(t — > /i7r°) arise from / ff and are very small. Although B(t — > /i?/, /i^'- 1 ) involve the 
poorly known parameters /i* and /iL , the two decays are dominated by ft^ and h^, , respectively. So the hadronic 
uncertainties are moderate. The upper horizonal lines denote current Belle upper limits. As shown in Fig.l, the 
branching ratios, as functions of |A TM |, rise rapidly with the increase of |A TAt |. We find that t~ — * fj,~rj(r]') decays are 

n 

more sensitive to |A TjU | than r — > /i ir° decay. From the Belle measurement |4j , we get bounds on the LFV coupling 
\X Tfi \ which are listed in Tabled for hadronic inputs with central values defaulted. t~ — > ii~r\,r<( decays give upper 
bounds on the strength of \X TiL \ at order of 0(1), while t~ — > [i~ir gives a looser bound. As shown in Tabled, 
comparing to the former constraints from other processes, we have obtained much more stringent bound for |A T/i | 
from r — » /i^''. Through calculations, we find that the branching ratios of t~ — > /i~7r can be as low as 2.0 x 10~ 9 
when X Tfl ~ 4.5. It should be noted that the current Belle bounds are based on only 153. 8/6 -1 data. Up to now, 
Belle has accumulated about 540/6~ 1 data already. To show the potential of bounding on LFV couplings with the 
full data at Belle, we scale the currents upper bounds naively by a factor 153.8/6 _1 /540/6 _1 . As a benchmark, the 
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potential are presented by lower horizontal lines in Fig.l, which would give bounds few times more restrictive than 
these from current Belle upper limits. 

To conclude this subsection, we have shown that, at the similar experimental sensitivity to the three decay modes, 
searching for r — > nrf would put more tighten constraints on the Higgs couplings than these from r — > (in decays. 



t -> fiM in RPV SUSY model 



The Supersymmetry model with explicit R-Parity breaking provides a simple framework for neutrino masses and 



mixing angles in agreement with the experimental data 



24| . The R parity quantum number is defined by 



3 



R = (-l) 



3B+L+2S 



(26) 



where B is the baryon number, L the lepton number and S the spin, respectively. 

Apparently, the lepton and/or baryon number violation could lead to R-Parity violation. The explicit R-Parity 
breaking would introduce renormalizablc bilinear higgsino-lcpton field mixings and trilinear Yukawa couplings between 
fermions and their super-partners 



(27) 



Lrpv - ^2 Hi{viR,HuL ~ e iRHuD + ^ { ^ijk 



A; 



1." 

2 



ViLdkRdjL + djLdkRViL + d^R^ijidjL — tihdkRUjL — UjLdkR&iL — d\ R e\ R Ujh 



" o\jk e af3~{ 



iaRd]PRdk 7 L + djpRUiaRd% lL + dl lR U iaR dj lL 



dL 



(28) 



where the indices i,j 1 k(= 1,2,3) label quark and lepton generations. Li and Qi are the SU(2)-doublet lepton and 
quark superfields. U^D^Ef are the singlet superfields. Ay* is antisymmetric under the interchange of the first two 
SU(2) indices, while X i - k is antisymmetric under the interchange of the last two. For the detailed reviews on R-parity 
violation SUSY model could be found in 3^, 3tl 4tJ 41 1. The phenomenological constraints on R- par ity violation 
couplings have been studied extensively in various processes with different points of view |42j, l43l |4J, |43, 1461 147 1 . For 
comparison, we list the constraints on A, A' couplings relevant to r — > fj,M in Table II, where the upper limits on A 
and A' are generally O(10~ 2 ) expect constraints from dd — ► uue~e~ and K — > ttvv. 
In the RPV SUSY model, the amplitudes of these decays are calculated to be 



A(t- -> h-tt°) = 
A(r~ -> fjL~rf) = 



1 



8V2'" , I 2m d 

1 (tfft 



'-L 2 (ht) 



Pi] f 7] 



V + A 



(iit) V -a 



h q 
"1; 



2(m u + m d ) 



2m 



11 T 



(fir) 



S-P 



(29) 



(30) 
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FIG. 2: B(t — > fin) as function of R-parity violation couplings, (a) for | A^ lx A*23 1 5 (b) for | A^ifc A 3 i fe | and | A^j! A 2 *i | ■ Others are 
the same as in Fig.l 



with 



A(T ->• H Tf) = - 



A' ;11 A.( 2 3 



1 \ P ' f ' 



Li Out) 



La 



P<J„' L ± o — L 3 



rj '7] 



(fir) 



V-A 



h q , 



2(m u + m d ) 



A 21fc A 31fc 



La 



h s , 

L\ + -^L 5 

2m, 



A 3jl A 2jl 



S-P 



L 4 



A 3j2 A 2j2 



A i22 A »23 



(31) 



(32) 



From Ea. (|29l31|l . we can see that the couplings relevant to lepton- number violation are A and A' couplings. There- 
fore, in the decays r — > fiM the parameters are the products of them: A 21fc A 3 i fc , AJ 11 Aj23, A^A^, ^j2^2j2^ an d 
A^ 22 Ai23 which are denoted by Li(i = 1,2,3,4,5), respectively. Among these coupling products, the first three 



TABLE II: Constraints on the products of A and A' R-parity violating couplings for sfermion masses ra± = lOOGeV . 



Couplings 



fin 



Previous Bounds 



|A-uA i2 3| < 1-46 x 10" 3 

lAai&Agifcl < 1.80 x 10~ 3 < 1.93 x 10~ 3 < 4.69 x 10~ 3 

I A^i Aa*-i | < 1.80 x 10~ 3 < 1.93 x 10~ 3 < 4.69 x 10~ 3 
|A3 j2 A2* 2 | < 9.38 x 10 -4 < 1.53 x 10~ 3 

|Ai22Ai23| 



< 3.6 x 10 _2 [45J, < 1.1 x 10" 2 [4£] 
lAjuA^il < 1.2 x 10- 2 [46] 
X 3j A < 9.1 x 10- 2 f47] 
lA^A^I < 1.2 x lO" 2 [46] 

< 4.5 x 10" 2 [45], < 3.0 x 10" 2 [48] 
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FIG. 3: B(t — > fir]) as function of R-parity violation couplings, (a) for A^A^I, (b) for | A^ife A 3 ^ fe and [A3, 1 A 2 5i|, (c) for 
|A3 j2 A 2 *2|, an d (d) for | Ai 2 2^i23 1 , respectively. Others are the same as in Fig.l 

contribute to r — ► /i7r° decay and all of them contribute to r — > (J,r)(r)') decays. Besides the type of (pr)s-p oper- 
ators, (J1t)v±a operators also appear in the amplitudes of these decays due to Fierz re-arrangements. Because of 
(P°\q^^q\0) = 0, A2 lfc A 3 * lfc and A^A^ will contribute to these decays with the same coefficient. 

As in literature, we assume that only one sfermion contributes one time with universal mass lOOGeV. We present 
our results in Fig. 2-4 for r _ — > /i~7r°, t~ — > /i~r; and t~ — * pTr( decays, respectively. 

From Fig. 2-4, we find that these decays could be enhanced to the present Belle sensitivities with the presences of 
RPV couplings constrained by other processes thus we can derive tighter bounds on the relevant RPV couplings. 
In Table II, we present bounds on RPV couplings derived from the Belle upper limits on the three r LFV decay 
modes|4| at 90% CL with central values defaulted for hadronic inputs. Most of them are stronger than before. 

For t — > /^7T° decay, only three coupling products |A^ n Ai23|, l^2ifc-^3ifc I ano - I^J5ji-^2ji I contribute to branching ratios. 
As shown by Fig. 2, the upper bounds for |A^ 11 Aj23|, lA^Ag*,.! and | A^ ^ A 2 * ^ | are 0(1O~ 3 ) which are more stringent 
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FIG. 4: B(t — > fir]') as function of R-parity violation couplings, (a) for |AinA;23|, (b) for | A^ife A 3 i fe and lA^ji A 2 **i |, (c) for 
1^3^2-^2*2! 1 an d (d) for I A^22 A,23 1 , respectively. Others are the same as in Fig.l 

than previous bounds (O(10- 2 ))Q |47| . It is noted that the contribution of |A£ u At2 3 | is enhanced by a factor 1/ ma- 
in, numerical calculation, we take ma, = (4.2 ± 1.0)MeV ^^|, as shown by Fig. 2(a), which causes large uncertainties 
for our theoretical prediction of |AJ 11 Ai23| contribution to t — > fin . With defaulted value = 4.2MeV, we get 
|^ u A i23 | < 1.46 x 10- 3 . 

All the five PRV coupling products contributing tor-> firf 1 fir}' decays. The sensitivities of the decays r — > fiif and 
t — » firf to these five RPV coupling products are depicted by Fig. 3 and Fig. 4, respectively. As shown by Fig. 3. (a) 
and Fig. 4. (a), the contributions of |AJ 11 Ai23|(ii) are subjected to huge theoretical uncertainties which arise from the 
poorly known ,^ and m q as shown by Eq. (|30l 131(1 . We could not get meaningful upper bounds on |A^ n Ai23| from 
t — > fjLTf^ fLv( decays. 

As shown by Fig.3.b-d and Fig.4.b-d, r — > firj,firf' decays are very sensitive to the contributions of |A2 lfe A 3 * lfe |, 
I ^371 ^2ji I > I ^3 1 2 -A>2^'2 1 an( i I A'-? ? A,-*q I ■ T herefor we get strong bounds on these four products which have improved the 



existing ones by one order [4f 



461 1471 l47y . From Fig.3.b-d and Fig.4.b-d, we can see that theoretical uncertainties are 
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small for predicting the contributions of these four RPV coupling products. 

It is noted that our study of r — > fir)' in RPV SUSY is new. The decays r — > »7T°, un have been studied by Kim, 

n n 

Ko and Lee|45j, however, we have used up-to-date hadronic inputs for 77 and 77' 30] . 

D. r — > fxM decays in Z 1 model with family non-universal couplings 

Many extensions of the standard model, especially grand unified theories and supersymmetry models, have addi- 
tional Z 1 bosons. In models with an extra U(l)' gauge boson, the Z' family non-universal couplings with the SM 
fermions generally induce flavor-changing neutral currents. In this paper, we refer to the basic formalism of Z' model 
elaborated in Ref. 

In the gauge basis, the Z' neutral-current Lagrangian can be written as 

C z ' = -g'J'^, (33) 

with g'the gauge coupling constant of the 1/(1)' group at the Mw scale. Here the renormalization group running 
between the Mw and Mz> scales is neglected considering the uncertainties of parameters. We assume that the Z' 
boson has no mixing with the SM Z boson The Z' current can be written as 

J l=zZ IMv P l + Ma?*] <t>h ( 34 ) 

i,3 

where / denotes the gauge interaction eigenstates and £^ L tm refers to the left(right)-handed chiral coupling matrix. 
The fermion Yukawa matrices in the weak eigenstate basis can be diagonalized by unitary matrices V ( / >L R 

yf a9 = V, a Y,VJ L . (35) 

V<I>l r could transform (j) 1 into mass eigenstate fields <Pl,r = V$ L r Pl,r4' 1 ■ The CKM matrix is given by the combi- 
nation 

Vckm = V U , L VI L (36) 

Flavor changing effects(FCNC) will present when e^ t rjA are nondiagonal matrices. The chiral Z' coupling matrices 
in the physical basis of fermions thus read 

= V uX e uX V^ x , B* = V dX e dX vJ x , (X = L, R) (37) 

where B*^ are hermitian. 

The low energy effective Hamiltonian of t — > fiM decays medicated by Z' is 

Hf /f - [BU^ p Lr) + B%{^P R r)\Y,[BU^ P ^)]+ h ^ ( 38 ) 
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where g\ = e/(sm9\y cos&w)- The diagonal elements of B L,R are real for hermiticity of the effective Hamiltonian, 
but the off-diagonal elements may contain weak phases. We introduce new weak phases (p L,R for B^ R {B^' R — 



\B 



L,R\ e i<f> • ) un d er assumption of neglecting B R . Now we can write out the decay amplitudes 



T/l 



TABLE III: Constraints on the values of £ in the flavor changing Z' model with family nonuniversal couplings 



Parameters 



flTT 



fig 



fig 



Previous Bounds 



\£ LL \ < 1.81 x 10" 3 < 1.91 x 10 -3 < 6.91 x 10" 3 < 0.02([-70°, -55°]), < 0.005([-80°, -30°]) [50]. 

< 0.0055(110°), < 0.0098(-97°) [51] 



\e 



< 1.81 x 10" 3 < 1.91 x 10" 3 < 6.91 x 10~ 3 



\^=^\ < 1.28 x 10" 



<1.34xl0" 3 <4.78xl0~ 3 < 0.02([5°, 15°]), < 0.005([-80°, -30°]) [5J 



< 0.0104(-70°), < 0.0186(83°) [51] 



G F 

A(T~ -> H~TT°) = -i — f^ap^o {X(p,T) V -A + Y(flT) V+A } , 



A(t n rj) 
A(t- -» iTrf) 



.Gp 



with 



x = / g'M z 
\9iM z , 

g'M z 



Ai = 

r x = 



As in literature, we take 



g'M z 

giM z > 

3 l&l 



pL r>L r>L r>L 



tdL r>L I pL tjL 

TfA UU + B Ttl B dd I 1 



y = 



(g'M z 



\9iMz 



r>R nl r-R r>£ 1 
n Tu n uu ~ n Tu n dd I ' 



A, 



g'Mz \ rL rL 



pi? pL I pi? pL \ p 



fj, JD ss' 



and £?£, 



^ B dd{s 



g' M Z N 2 pL R L pR£ 



(39) 
(40) 
(41) 

(42) 
(43) 
(44) 

(45) 



We find that the variation of phase have nearly negligible influence on the branching ratios. So we take 

Pl,r = io° Q as benchmark. At first, we assume that only one parameter contributes one time and then both £ 
and £ RR contribute with the same strength. We present the correlation of branching ratio versus parameter |£| in 
Fig. 5-7. The branching ratios get the same contributions from \£ LL \ and |£ |. The upper bounds on parameters 
|£| from experimental data are 0(1CP 3 ) which are listed in Table III for comparison. Compared to former two new 
physics scenarios, there is no scalar operator induced by Z' family non-universal couplings. Thus the poorly known 
hadronic parameters h q ^ and factor \/m q are absent. Theoretical uncertainties are due to decay constants f v and 
/^'*/ which have been determined with few percentages accuracy as listed in Eq.(10). So theoretical estimations of 
these decays could be made quite accurate. 
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FIG. 5: B(t -> fin ) as function of Z' couplings, (a) for \£ LL \ and \£ RL \ contributing seperatly. (b) for both \£ LL \ and \£ RL \ 
contributing with the same strength. The horizontal lines are the same as these in Fig.l. 




FIG. 6: B(t — > fir]) as function of Z' couplings. The horizontal lines are same as these in Fig.l. 

As shown by Fig. 5 and 6, the sensitivities of B(t — ► jU7T°) and B(t — > /wr ) to the Z' LFV FCNC couplings are quite 
similar. From Belle upper limit on B(r — > \i% ), we get ( < 2x 10 -3 which are comparable with previous bounds 
[50! 51 1 . However, from Fig. 7, we can see that r — * fit]' decay gives weaker bounds than r — ► pkr)(iT). 

Anyhow the Belle searching for r — > /jM decays have already put strong bounds on the parameter spaces. In other 
words, these decays could be enhanced to the present B factories sensitives by the Z' family nonuniversal FCNC 
couplings without conflict with bounds from other exist measurements. 
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5x1(T 1(T 

Z* parameter 

(a) 





„LL „RL 










3x1 (T 

Z parameter 



6x1 <T 



FIG. 7: B(t — > /zr/) as function of Z 1 couplings. The horizontal lines are same as these in Fig.l. 



III. CONCLUSION 



The measurement of LFV processes would be a definite evidence for physics beyond the SM. In this paper, we 
have studied LFV processes t" -> fi~M(M = tt°, 77, 77') at the tree level in the 2HDM III, RPV SUSY and flavor 
changing Z 1 model with family non-universal couplings. Since these decays are very sensitive to the presences of LFV 
couplings, we have derived constraints on parameter space of the three New Physics scenarios from the recent Belle 
limits 4]. Our main findings can be summarized as follows. 

1. In 2HDM III, the strongest bound on X Tfi comes from r — > [ir\ decay. The bound is consistent with these in 
literature, however, improve these by several times. 

2. In RPV SUSY, r — * ^77,^7/ decays are very sensitive to the contributions of |A' 21fc A 3 * lfe |, (A^A^'il) I^3j2-^2j"2 1 
and |A^ 22 Aj23l Therefor we g et strong bounds on these three products which have improved the existing ones 
by one order [ifl 4fl 42, 47J. However, there are large uncertainties in calculating contributions of these RPV 
coupling products due to hadronic parameters ,. We could not get bounds on lA^A^I from the two decays 
because of poorly known h s v _/. However, we can get strong bound on lA^A^gl with r — * /U7r°. 

3. In Z' model, theoretical predictions of Z' LFV coupling contributions could be make quite accurate, r — > /J7r, fir] 
and fir] 1 decays have similar sensitivities to the LFV couplings. Belle current upper limits on B(r — > fj,i], fir]') 
have already constrained £ as small as 0(1CP 3 ). 

In summary, we have shown that the LFV semilcptonic decays r — > fin , /irj and /i?/ are very sensitive to the 
presences of LFV couplings in 2HDM III, RPV SUSY and Z' models. Using only 153. 8/6 -1 data, Belle recent upper 
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limits for these decays have already given quite tight bounds on the LFV couplings in the aforementioned three new 
physics scenarios. It should be noted that Belle and BaBar had accumulated about 540/6 -1 and 320/6 -1 , respectively, 
till the end of year 2005. With refined measurements with the these data at Belle and BaBar, we could get more 
crucial information on LFV, at least more stringent bounds on the parameter spaces of models with LFV couplings. 
The results derived in this paper from the recent measurements at Belle would be useful for phenomenological studies 
of the scenarios in other interesting processes. 
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